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Abstract 

In this work a simple problem on 2D optimal shape for body im- 
mersed in a viscous flow is analyzed. The body has geometrical con- 
straints and its profile would be found in the class of cubics which 
satisfy those conditions. The optimal profile depends on the leading 
coefficient of these cubics and its relation with the Reynolds number 
of the system is found. The solution to the problem uses a method 
based on a suitable transformation rule for the cartesian reference. 
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We would study the problem of searching the optimal 2D shape or profile 
of a cartesian object immersed in a constant fluid viscous flow and subjected 
to some constraints on the boundary. This is a particular question on the 
general field of shape optimization for bodies immersed in flows (see [1]). 

Let be {x,y} a cartesian system, and y = f{x) a function whose graph 
is the profile which we want to optimize. The function / is subjected to the 
following constraints, arising from engineering requirements: 
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cl. /(O) = 

c2. /(xo) = yo 

c3. f'ixo) = 

c4. f(0) > 

c5. fix) > Vx G (0,xo) 

where Xq and ?/o ^^re positive numbers. The graph is immersed in a vis- 
cous flow having, at freestream zone, constant velocity V = {voo,0) parallel 
to the X-axis, with Voo > 0. What is, or which are, the functions / that 
minimize the pressure drop along its profile, that is, if p = p{x, y) is the 
static pressure along the graph, what is the curve for which the difference 
Ap = p(0, 0) — p(xo, yo) is minimum? The question is related to drag and lift 
optimization (see [T]). 

We consider the cubics /(x) = ax^ + 6x^ + cx + d. Using cl, c2 and c3, 
we find 



yo + 2axl 



Xq 

22/0 + axl . 
c = (1) 

Xo 

d = 

The derivative of the general cubic is /'(x) = 3ax^ + 2fex + c so that /'(O) = c. 
From c4, we have the following condition on parameter a: 

a > -2^ (2) 

while a few algebra, imposing the condition f'{t) = if and only if t < or 
t > Xq, shows that c5 is satisfied if 

^ 2/0 



Then the search of the optimal profile is defined on the class of cubics 



^^^^^,_y_o+pl^,^2yQ + a4^ . < , < ^ j (4) 

Xq xq Xq Xq ) 



2 



y 




0.5 1.0 1.5 2.0 



Figure 1: Some cubics in the case Xq = 2, uq = 1: from a = —0.25 (blue 
curve) to a = 0.125 (green). 



We use Newton's theory on fluid velocity distribution on the profile or 
surface of a body immersed in a flow with constant freestream speed 
([2]). If a = a (x) is the angle between the direction of the flow and the 
tangent to the curve, the velocity in a point of the profile has two components, 
the tangential one with module VooCosa, and the normal one with module 
VooSina. The latter gives the amount for the body resistence (see [3]). But 
we would analyze the flow on the upper neighbour of the profile, where fluid 
has a tangent velocity field given by 

V = (voocos^a, VooSina cosa) (5) 
From usual trigonometrical formulas the following identity holds 

cosM^) = (6) 

so that we have 

Let fi the dynamic viscosity, p the density and p = p{x, f{x)) the static 
pressure of the fluid on the body profile. We could write the Navier-Stokes 
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equations for this flow in the upper neighbour of the profile, with ([7j) and 
p{0, 0) = p{xo, Uo) as boundary condition: 



p Vv V = —Vp + /iAv 

p(0,0) =p{xo,yo) 
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) 



We try to simplify the resolution of previous system by the following trans- 
formation rule on coordinates system: 



It is important to note that, from c5 condition and usual notions of real 
analysis, the function / is invertible on the interval (0,a;o) with inverse g 



In the new coordinates system {X,Y}, the curve y = f{x) has the simple 
equation Y = 0, hence it becomes the X-axis itself. 



Figure 2: The new cartesian system XY , from the old xy point of view, in 
the case a = —0.25. If X = and y = 0, then x = and Y = so that the 
new origin is the same as the old. For Y = we have y = f{x), so that the 
new X-axis is the curve f in the old system. If X = 0, then f{x) = and 
Y = y, so that the new Y-axis is the old y-axis in the xy system. The thin 
line has equation of type Y = cost. 

Then, the (X,y)-fiow near the curve is parallel to X, that is the velocity 
field U in the system {X, F } has only the X-component: U = {Ui,0). If 
a = p{giX),Y + X),u = fiigiX),Y + X) and P = p{g{X), Y + X) are the 




(8) 



differentiable: dxg = j-j- 



y 




X 
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representation of the scalar functions p, /i, p in the {X, Y} system, the flow 
is described by the simple Navier-Stokes equation 

dU, dP d'U, 
The boundary conditions are now U{X, 0) = vi{g{X), Y + X) and P(0, 0) = 

Piyo,o). 

We can write the expression of Ui, 

U, = U,iX,Y) = vMX),Y + X) = VooY^ (10) 

where g' = dxg, and the expression of its first derivative: 

mh_ 2g'g"{l + g'^)-2g'^g" _^ g'g" 

dX {l + g'^y + ^^^^ 

Note that Ui doesn't depend explicitly on Y. For X = /(xq) = yo and Y = 0, 
the value of Ui is defined by continuity extension, because /'(xq) = and 
consequently dxgivo) = oo- From (flOl) 

c?'2 1 

Uiivo, 0) = lim ^;oo— — ^ = ^^oo hm = v^o (12) 

A^yo I + g x^yo ^ + i 

In the case /'(O) = 0, at the same manner the definition of Ui can be extended 
at (0,0) by U,{0,0) = v^. 

Now we integrate the two members of respect to the X variable from 
to f{xo) = yo: 

l^[un7 = -[P]o+^[dxUi]l° (13) 

For our purpose of optimization, we impose — [P]q° = F(0,0) — P{yQ,0) = 
p(0,0) -pixo,yo) = 0. 

At first consider the simple case of inviscid fiow: u = 0. Previous equation 
becomes 

l-[un:=^ (14) 

therefore it must be ?7i(0,0) = Ui{yQ,0). Using previous identities, this 
equation can be written as 
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(15) 



(i + /'(o)T. 

Then it must be /'(O) = 0, that is c = = 0. Therefore, in the case of 

inviscid flow, the optimization problem is solved by the cubic, which belongs 
to the class C, with 



-2^ 



(16) 



The cartesian equation is y 





-2^x3 + 3^x2. 



Now we consider the viscous case. The condition to impose in equation 
(IT^ is always [P]q° = 0, therefore 



(17) 



The first member is 



For expliciting the second member, from (fTTl) we have to find g" = gxx- 
Apply usual differentiation rules: 



gxxiX) = dxgx{X) = d 
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X 



1 /" 

d^cjdxx = -j^gx 



dJ{g{X)) 

f"{g{X))gUX) 



(19) 



Therefore we can write 
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dxUi{yo) = - Jim 2v^f\g{X)) ^ = -2v^f"{xo) (20) 



dxUM = - lim^2v^ng{X))j^^i2L_^ = -2v^ ,''' (21) 



xo ' 

/"(O) = _2^2±^ and /"(xq) = 2 '^°^""° ■ Equation can be written as 



For a cubic of the class C the following identities hold: /'(O) 

-yo+ax^ 



avr, 



no) 



(22) 



or 



av^ [{1 + f'iOff -l]=Au [f"{0) - (1 + fm^fixo)] (23) 

It is an algebraic equation of the 5th order in the parameter a. But we are 
interested in physical situations where viscosity is small, that is when the 
parameter a has values near the previous computed quantity fll6p : then, in 
this situation, /'(O) is small, /"(O) = 2fe = 6|| and f"{xo) = -6g = -/"(O). 
We can consider the following Taylor expansion of the quantity (1 + /'(O)^)^: 



;i + .no) 



2\2 



1 + 2/'2(0) = 1 + 2 



{2yo + axl)^ 



Xn 



(24) 



Then the algebraic equation in the parameter a can be simplified in the 
following 2nd order one: 



1 + 2^^y^±^ ] {Auf"{0) - av^) + (4z./"(0) + av^) = (25) 

Xr\ 



The admissible solution to this equations is 



-2^ + 1 



Xn 



Xn 



|24z/?/o - crvoox, 



(26) 



Note that, in the case = 0, the solution has the expression (1161) . as expected. 
In viscous case, on the contrary as inviscid one, the optimal profile depends 
on the values of density p, viscosity p, and freestream speed Voo- 
We can see that the main effect of viscosity is its influence on /'(O), that is 
on the angle of attack between flow and profile. In fact, a is an increasing 



function of p, therefore /'(O) = c 



2yo+ax^ 



xo 



is an increasing function too. 



If we consider a flow with a constant value of /i, it is interesting to evaluate 
a as function of the freestream speed Voo- From fl26l) follows that a is a 
decreasing function of f oo and 



lim a{vc 



Xn 



(27) 
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aOu) 



-0.248 



-0.242 



-0.244 



-0.246 
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Figure 3: Parameter a as function of jjL in the case xq = 2, uq = 1, a = 0.001 
and Voo = 1000 {CGS units). 



therefore the increasing of flow speed is equivalent to a vanishing of viscosity. 

Now multiply by j/o and divide by vx^ both numerator and denominator 
of the radicand in fl26l) . Introducing the label = |^ and the Reynolds 
number 



(recall that a is the density and v the dynamic viscosity) where yo is the 
characteristic length of this geometrical system (see [1]), the parameter a 
can be written in the form 



This expression separates the dependence of a on geometrical (xq, 2^0) and 
physical parameters (Me). As expected, for high Reynolds number the influ- 
ence of viscosity vanishes and the optimal profile tends to the shape of the 
case /i = 0: 




(29) 



Re->oo X\ 



lim a(Re) = —{-2zo) 




(30) 
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a([Re) 



-0.000058 
-0.0000585 ; 

-0.000059 
-0.0000595 



Re 

20000 40000 60000 80000 100000 

Figure 4: Parameter a as function of Re in the case Xq — 100, Zq — 0.3 
{CGS units). 
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Figure 5: Optimal profile as function of /j, & [0.0, 0.001] in the case xq — 
100; zq — 0.3 {CGS units); decreasing grey color levels are associated to 
increasing values of viscosity. 
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